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Abstract 

The method of the quantum probability theory only requires simple structural data of 
graph and allows us to avoid a heavy combinational argument often necessary to obtain 
full description of spectrum of the adjacency matrix. In the present paper, by using the 
idea of calculation of the probability amplitudes for continuous-time quantum walk in 
terms of the quantum probability theory, we investigate quantum central limit theorem 
for continuous-time quantum walks on odd graphs. 
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Two types of quantum walks, discrete and continuous time, were introduced as the quantum 
mechanical extension of the corresponding random walks and have been extensively studied 
over the last few years [1, 2]. 

Random walks on graphs are the basis of a number of classical algorithms. Examples include 
2-SAT (satisfiability for certain types of Boolean formulas), graph connectivity, and finding 
satisfying assignments for Boolean formulas. It is this success of random walks that motivated 
the study of their quantum analogs in order to explore whether they might extend the set of 
quantum algorithms. Several systems have been proposed as candidates to implement quantum 
random walks. These proposals include atoms trapped in optical lattices [3], cavity quantum 
electrodynamics (CQED) [4] and nuclear magnetic resonance (NMR) in solid substrates [5, 6]. 
In hquid-state NMR systems [7] , timeresolved observations of spin waves has been done [8] . It 
has also been pointed out that a quantum walk can be simulated using classical waves instead 
of matter waves [9, 10]. 

A study of quantum walks on simple graph is well known in physics(see [11]). Recent 
studies of quantum walks on more general graphs were described in [12, 13, 14, 15, 16, 17, 18]. 
Some of these works studies the problem in the important context of algorithmic problems on 
graphs and suggests that quantum walks is a promising algorithmic technique for designing 
future quantum algorithms. 

There is one approach for investigation of continuous-time quantum walk (CTQW) on 
graphs using the spectral distribution associated with the adjacency matrix of graphs [19, 20, 
21, 22, 23]. Authors in Ref.[19] have introduced a new method for calculating the probability 
amplitudes of quantum walk based on spectral distribution which allows us to avoid a heavy 
combinational argument often necessary to obtain full description of spectrum of the Hamilto- 
nian. It is interesting to investigate the CTQW on graph when the graph grows as time goes 
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by. In order to study a quantum system in full detail, its Hamiltonian needs to be diagonalized. 
With increasing dimension of the Hilbert space, the diagonalization of an operator becomes a 
very tedious task. In fact, we discuss this question in the CTQW as a quantum central limit 
theorem for CTQW. In this paper we try to investigate quantum central limit theorem for 
CTQW on growing odd graph via spectral distribution. 

The organization of the paper is as follows. In section 2, we give a brief outline of graphs and 
introduce odd graph. In Section 3, we review the stratification and quantum decomposition for 
adjacency matrix of graphs. Section 4 is devoted to the method of computing the amplitude 
for CTQW, through spectral distribution of the adjacency matrix. In the first subsection 
of section 5 we evaluate the CTQW on finite odd graph and the in second subsection 5, we 
investigate quantum central limit theorem for CTQW on odd graph. The paper is ended with 
a brief conclusion and one appendix which contains determination of spectral distribution by 
continued fractions method. 



In beginning we present a summery of graph and then we introduce odd graph. 

A graph is a pair G = (y,E), where F is a non-empty set and is a subset oi {{a, P}\a, (3 e 
V,a 7^ /?}. Elements of V and E are called vertices and edges, respectively. Two vertices 
a,P eV are called adjacent if {a, P} e E, and in this case we write a ~ /3. Let l'^{V) denotes 
the Hilbert space of C-valued squarc-summable functions on V, and {\a)\ a G V} becomes a 
complete orthonormal basis of l'^{V). The adjacency matrix A — {Aai3)a,i3ev is defined by 



2 Odd graph 



1 if a ~ 




otherwise. 



which is considered as an operator acting on P 



(V) in such a way that 



= E \P)^ 



aeV. 
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Obviously, (i) A is symmetric; (ii) an element of A takes a value in {0, 1}; (iii) a diagonal 
element of A vanishes. Conversely, for a non-empty set V, a graph structure is uniquely 
determined by such a matrix indexed by V. The degree or valency of a vertex a &V is defined 

by 

K(a) ^ \{(3 e V\a (3}\, 

where |.| denotes the cardinality. 

Let 5" be a set of integer as 5" = {1, 2, 2k — 1} for a fixed integer k > 2. Now we define 
V be the set of subsets of S having cardinality k — 1, i.e., 

V ^{aG S \a\^k-l}, (2-1) 

and put 

E = {{a,P}\a,P eV,anp = (l)}. (2-2) 

This graph (V, E) is called the odd graph of degree k and is denoted by Ok- Obviously, Ok is a 
regular graph with degree k. For k = 3, the odd graph is well-known as Petersen graph (Fig.l). 
The odd graphs have been studied in algebraic graph theory where some of their properties 
are found in Refs.[24, 25]. 

For n = 0, 1, 2, A; — 1 we define e„ as 

A; — 1 — I if n is even 
^ if n is odd, 

then by using proposition 4.1 [26], for pair a, P E V we have 

|Q;n/3| = e„ <^ 9(q;, /3) = n, (2-3) 

where d stands for the natural distance function. Due to the above relations based on distance 
function the odd graphs are distance regular graphs, i.e., for a given — 0,1,2,..., the 
intersection number 

p^,. = |{7ey|5(a,7)=^and d{^,P)=j}\, (2-4) 
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is independently determined of the choice of Q!,/3 e V satisfying d{a,(3) = I. There are some 
well-known facts about the intersection numbers of distance regular graphs, for example ^ = 
(for i 7^ 0, j is not {i — + 1})( for more details see Refs. [27, 28]). For convenience, set 
bi := <^<d),Ci:= pl^,^,{0 < z < (i - 1), a, := pl,{0 <i<d),h:= pO,(0 < i < d) 

and hQ — Cd — 0, where d :=max{9(Q;, (3) : a, (3 is called diameter of graph. Moreover 

hi + ai + Ci^k, 0<i <d, (2-5) 

where k — ki is degree of graph. 

For calculation of Szego- Jacobi sequences {0;^} and {ctfc} in the next section we will require 
bi and q [20], where for odd graph is given by 

{I if i is even, 
(1<^<A;-1), (2-6) 
^ if i is odd 

{/c — I if i is even, 
{0<i<k-2). (2-7) 
A; - i±i if i is odd 

3 Quantum Probabilistic Approach for CTQW 

In this section we give some preliminaries that require to describe CTQW via spectral distri- 
bution technique. 

3.1 Stratification 

Due to definition of function 9, the graph becomes a metric space with the distance d. We fix 
a point o e y as an origin of the graph, called reference vertex. Then, the graph is stratified 
into a disjoint union of strata: 

00 

V^lJVi, Vi^{ae V\ d{o, a) = i}. (3-8) 

1=0 
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With each stratum Vi we associate a unit vector in l'^{V) defined by 

1 



|0i 



(3-9) 



\Vi\ aeVi 

where \i,a) denotes the eigenket of the ath vertex at the stratum i. The closed subspace of 
l'^{V) spanned by is denoted by A(G). Since becomes a complete orthonormal 

basis of A(G), we often write 

A(G)=^©C|0,). (3-10) 



3.2 Quantum decomposition 

One can obtain a quantum decomposition associated with the stratification (3-8) for the ad- 
jacency matrices of this type of graphs as 



A^A+ + A- + A°. 



where three matrices ^4+, A and ^4° are defined as follows: for o; e 



if /3 e Vi^x 

otherwise, 



(^")/3a = < 



(^°)/3a = < 



An^ if/Sey^-i 

otherwise, 

A^^ if/Se^i 

otherwise. 



or, equivalently, for |i, oi) 



(3-11) 



N + A-\i,a)^ a) = ^ |i, /?), (3-12) 

for {a,/?} e -B. Since a e Vi and {a,/?} e E then /3 e V^-i U U ^i+i, where we tacitly 
understand that V^i — 0. The vector state corresponding to |o) = |0o)) with a e V as the 
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fixed origin, is analogous to the vacuum state in Fock space. According to Ref.[29], < A"^ > 
coincides with the number of m-step walks starting and terminating at a, also, by lemma 2.2, 
[29] if G is invariant under the quantum components A'^, e G {+,—,0}, then there exist two 
Szego- Jacobi sequences {uji}'?!^ and {0;^}^]^ derived from A, such that 

= I > (3-13) 

A-\(Po)^0, = v^|0,_i), i>l (3-14) 

= ^>0, (3-15) 

where ^^uJi+i = '^^ii/T ^-(/S), = |{a e Vi\a ~ for /? G and a^+i = ko(/3), such 

that «;o(/3) = £ ~ /?}| for G V^. In particular {A{G), A'^ , A~) is an interacting Fock 
space associated with a Szego- Jacobi sequence {cui}. 

3.3 Study of CTQW on a graph via spectral distribution of its 
adjacency matrix 

The CTQW on graph has been introduced as the quantum mechanical analogue of its clas- 
sical counterpart, which is defined by replacing Kolmogorov's equation (master equation) of 
continuous-time classical random walk on a graph [30, 31] 

^%^ = ^i/„„P,(i), m=l,2,...,/ (3-16) 

n=l 

with Schrodinger's equation. Matrix H is the Hamiltonian of the walk and Pm{t) is the 
occupying probability of vertex m at time t. It is natural to choose the Laplacian of the graph, 
defined as L — A — D as the Hamiltonian of the walk, where D is a diagonal matrix with 
entries Djj = deg{aj). 

Let \(p{t)) be a time-dependent amplitude of the quantum process on graph F. The wave 
evolution of the quantum walk is 

ih^^\m-H\m), (3-17) 
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where from now on we assume h = 1, and \4>q) is the initial amphtude wave function of the 
particle. The solution is given by \(f){t)) — e~*^*|0o)- On d-regular graphs, D — |7, and since 
A and D commute, we get 

This introduces an irrelevant phase factor in the wave evolution. Hence we can consider H — A. 
Thus, we have 

|0(t))=e-^^*|0o), (3-19) 

One of our goals in this paper is the evaluation of probability amplitudes for CTQW on 
graphs by using the method of spectral distribution associated with the adjacency matrix. The 
spectral properties of the adjacency matrix of a graph play an important role in many branches 
of mathematics and physics. The spectral distribution can be generalized in various ways. In 
this work, following Refs.[19, 29], we consider the spectral distribution ^ of the adjacency 
matrix A: 

iA"^) = / x'^iiidx), m = 0, 1, 2, ... (3-20) 
Jr 

where (.) is the mean value with respect to the state |0o)- By condition of quantum decomposi- 
tion (QD) graphs the moment sequence {(^"*)}m=o well-defined [19, 29]. Then the existence 
of a spectral distribution satisfying (3-20) is a consequence of Hamburgers theorem, see e.g., 
Shohat and Tamarkin [[33], Theorem 1.2]. 

Due to using the quantum decomposition relations (3-13, 3-14,3-15) and the recursion 
relation (A-i) of polynomial Pn{x), the other matrix elements {4'n\A™' \ 4>q) can be written as 

(</.„|^'" I 0o) = , ^ I x"'Pn{x)i^{dx), m = 0, 1, 2, .... (3-21) 

where is useful for obtaining of amplitudes of CTQW in terms of spectral distribution associ- 
ated with the adjacency matrix of graphs [19]. 
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Therefore by using (3-21), the probabihty amphtude of observing the walk at stratum m 
at time t can be obtained as 

qm{t) = = (0^|e-^^* I 0o) = , ^ / e-'^'Pm{x)pi{dx). (3-22) 

The conservation of probabihty Y,m.=G I (0m|0(^)) 1^ = 1 foUows immediately from Eq.(3-22) 
by using the completeness relation of orthogonal polynomials Pn{x). In the appendix A of the 
reference [19], it is proved that the walker has the same amplitude at the vertices belonging 
to the same stratum, i.e., we have qim{t) = = 0, 1, | |, where qim{t) denotes the 

amplitude of the walker at ith vertex of mth stratum. 

Investigation of CTQW via spectral distribution method pave the way to calculate CTQW 
on infinite graphs and to approximate with finite graphs and vice versa, simply via Gauss 
quadrature formula, where in cases of infinite graphs, one can study asymptotic behavior of 
walk at large enough times by using the method of stationary phase approximation (for more 
details see [1]). 

Indeed, the determination of jj,{x) is the main problem in the spectral theory of operators, 
where this is quite possible by using the continued fractions method, as it is explained in 
appendix A. 



4 Quantum central limit theorem for CTQW on odd 
graphs 

Having studied CTQW on finite odd graphs using the method of the spectral distribution, we 
investigate quantum central limit theorem for CTQW on this graphs which is our main goals. 

To consider stratification and quantum decomposition of section 3, Ref. [20] (i.e., cui = 
Ci-ibi, ai = ai — bi-i — Q-i) and Eqs.(2-6), (2-7), we obtain two Szego- Jacobi sequences 
{cUi} and {ccj} for as follows: 
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if i is odd, 



if i is even, 



u^^-'^ik-'-^), (4-23) 



^. = ^(^-^) (4-24) 



if i is odd or even, 

ai = 0. (4-25) 

A. Finite k case 

Let Ilk denote the spectral distribution of odd graph Ok- Here for studying CTQW on 
finite odd graph we consider k = 4 case. Then we have 

uji = i, a;2 = 3, LUs = 6, ai = 02 = • • • = 0. (4-26) 

Therefore we obtain Stieltjes transform 

In this case one can obtain the spectral distribution as follows 
II4 



292 

73, 



73 1 / 1 / 

(5 + V73){d{x - -V26 - 2^73) + 6{x + -\/26 - 2^73))+ 



5 + V73){S{x - V26 + 2V73) + S{x + V26 + 2^73)) (4-28) 



292 ' - /V V 2 * ' ' 2 

By using Eq.(3-22) the amplitudes for walk at time t are 



/73 1 / 1 / 

g„(i) = ^((5 + V73) cos(-V26 - 2V73)t + (-5 + ^73) cos(- V26 + 2^73)^), 

gi(t) = ^^^((5 + V73)(V26-2V73) sin(i\/26 - 2^73)^+ 

,1 



'73 



-5 + \/73)(V26 + 2\/73) sin(-V26 + 2V73)t), 
^(- cos(--\/26-2/73)t + cos(--\/26 + 2\/73)t). 
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qzit) = ^^^(-(13 + V73)(\/26 - 2V73) sm(^\/26 - 2^73)^+ 

(-13 + V73)(V26 + 2V73) sin( ^26 + 2V73)t). (4-29) 

B. Quantum central limit theorem 

In the limit of large k — > oo, it is observed the odd graphs form a growing family of 
distance regular graphs. In the remaining of this section, we obtain CTQW on odd graphs 
Ofc for k — > oo by applying the quantum central limit theorem where is derived from the 
quantum probabilistic techniques. 

Theorem ([26]). Let = {V'',E'^)} be a growing family of distance regular graphs. 
Let Ak and {Pij{k)} be the adjacency matrix and the intersection numbers of G^, respectively. 
Assume that the limits 

Pl-i,i{k)pl:^{k) h{k)c,_,{k) 
uji = hm n n\ = 1™ ; 

k — >oo Pll\k) — >oo k 

= Hm = lim ' h-^{k) - c^.,{k) 

exist for alH = 1, 2, .... Let A = (A, {[■^j)}, B~) be the interacting Fock space associated 
with {uJi} and define a diagonal operator by B^\ipi) = Q;j+i|?/'j). Therefor the quantum 
components A^., e e {+, — , 0}, of adjacency matrix A^ is holds that 

lim = lim 4^ = 5^ £ e {+, -, 0}, (4-31) 

y'p?i(A;) ''-^'^ Vk 

in the stochastic sense. Then we have 

lim (0^|4|l0o) = (^ml^'l^o), (4-32) 

for e e {+, -,0}. 

To state a quantum central limit theorem for CTQW on odd graph Ok, it is convenient to 
calculate amplitudes of probability as 

qm{t) = lim (0^|e^|0o) = , ^ / e-'^' Prn{x) fi^{dx) . (4-33) 

k — >oo ~^y0iJiLV2---'^m JR 
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According to the above theorem, we need only to find two Szego- Jacobi sequences {uui} and 
{ctj}, where by using Eqs.(4-30) and (4-23) we obtain as follows: 
if i is odd, 



if i is even. 



li + 1,, i — Is i + 1 /..X 



if i is odd or even, 

ai = 0. (4-36) 

Thus, {ui} = {1,1,2,2,3,3,4,4,...} as desired. Therefore Stieltjes transform of infinite odd 
graphs is 

G,^ (z) = ^ 1 (4-37) 



~5~ 

z 

where the spectral distribution ^oo{x) in the Stieltjes transform is given by [26] 

Atoo(x) = |x|e-^'. (4-38) 

Therefore we obtain the amplitude of probability at time t and the 0-th stratum (starting 
vertex) 

JR J-oo 

~ ~^ 

In the above calculation we used formulas: 



erf{it/2)e-^ + 1. (4-39) 



/ e"^ e-*^*dx= ^(l-er/(it/2))e— , (4-40) 

JQ 2 

where erf{x) stands for the error function and it is defined as: 

erf{x) = ^ re-'^ds (4-41) 



Jo 
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also, the derivative of the error function follows immediately from its definition ^erf{x) — 

By using Eqs.(3-22) and (4-40), we obtain the amplitude of probability for walk at time t 
and m-th stratum on infinite odd graphs in terms of qo{t) as: 
if m is odd 

Qm{t) = J^^T^^^Prn{tf^)qoit), (4-42) 

if m is even 

Qm{t) = (4-43) 

where polynomials {Pm{i-^)} are defined recurrently by relation (A-i) in terms of ij^. In the 
end for example we obtain the amplitude of probability for walk at stratum 1, 2 and 3 as 
follows: 
m — 1 



m = 2 



q,{t) = Pl{^f^)qo{t) = if^qo{t) = ^{t' - 2)erf{zt/2)e^ - zt/2, (4-44) 



q2{t) = ^2(^|)go(^) = - l)qo{t) 

^ ' iy/nt^ erf{it/2)e^ - 2t^ + 2iy/^t erf{it/2)e^), (4-45) 



m — 3 
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5 Conclusion 
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\/7rt^er/(it/2)e^ + 2it^ + A^t^erf{it/2)e'^ - Ait + A^/Tierf{it/2)e'^). (4-46) 



In this paper by using the method of calculation of the probability amplitude for continuous- 
time quantum walk on graph via quantum probability theory, we have studied continuous-time 
quantum walk on odd graph when the odd graph grow as time goes by. We have discussed 
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this question as a quantum central limit theorem for CTQW. It is interesting to investigate 
CTQW on a growing family of graphs since it is probable the probabihty amplitudes of CTQW 
to converge the uniform distribution, which is under investigation. 



Appendix A 

Determination of spectral distribution by continued fractions method 

In this appendix we explain how we can determine spectral distribution ii{x) of the graphs, 
by using the Szego-Jacobi sequences {{ujk}, {o^fc}), which the parameters Uk and ak are defined 
in the subsection 3.2. 

To this aim we may apply the canonical isomorphism from the interacting Fock space onto 
the closed linear span of the orthogonal polynomials determined by the Szego-Jacobi sequences 
{{uJi}, {o(i}). More precisely, the spectral distribution /i under question is characterized by the 
property of orthogonalizing the polynomials {Pn} defined recurrently by 

Po{x) = 1, Pi{x) = X - «!, 
xPn{x) = Pn+l{x) + an+lPn{x) + UJnPn~l{x) , (A-i) 

for n> 1. 

As it is shown in [32] , the spectral distribution can be determined by the following identity: 

G (z)= f ^^^^^ = - = ^n-ljz) ^ A 

' Jrz-x z-a,- ^ -2^ - Pn{z) ^iZ-xi' ^ ' 

^-°3- ^-ri-... 

where G^{z) is called the Stieltjes transform and Ai is the coefficient in the Gauss quadrature 
formula corresponding to the roots xi of polynomial Pn{x) and where the polynomials {Qn^} 
are defined recurrently as 

Q^i^ix) ^ X - a2, 
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for n > 1. 



Now if G^{z) is known, then the spectral distribution fj, can be recovered from Gfj,{z) by 
means of the Stieltjes inversion formula: 



Substituting the right hand side of ( A-ii) in (A-iii) , the spectral distribution can be determined 
in terms oi xi,l — 1,2, the roots of the polynomial and Guass quadrature constant 

Ai, I = 1,2, ... as 



( for more details see Refs. [19, 20, 32, 33].) 
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Figure Captions 

Figure- 1: The Petersen graph. 



